Development of an improved divergence‐free‐condition compensated coupled framework to solve flow problems with time‐varying geometries

In this article a coupled version of the improved divergence‐free‐condition compensated method will be proposed to simulate time‐varying geometries by direct forcing immersed boundary method. The proposed method can be seen as a quasi‐multi‐moment framework due to the fact that the momentum equations are discretized by both cell‐centered and cell‐face velocity. For simulating time‐varying geometries, a semi‐implicit iterative method is proposed for calculating the direct forcing terms. Treatments for suppressing spurious force oscillations, calculating drag/lift forces, and evaluating velocity and pressure for freshly cells will also be addressed. In order to show the applicability and accuracy, analytical as well as benchmark problems will be investigated by the present framework and compared with other numerical and experimental results.     

Analysis of a new multispecies tumor growth model coupling 3D phase-fields with a 1D vascular network

In this work, we present and analyze a mathematical model for tumor growth incorporating ECM erosion, interstitial flow, and the effect of vascular flow and nutrient transport. The model is of phase-field or diffused-interface type in which multiple phases of cell species and other constituents are separated by smooth evolving interfaces. The model involves a mesoscale version of Darcy’s law to capture the flow mechanism in the tissue matrix. Modeling flow and transport processes in the vasculature supplying the healthy and cancerous tissue, one-dimensional (1D) equations are considered. Since the models governing the transport and flow processes are defined together with cell species models on a three-dimensional (3D) domain, we obtain a 3D–1D coupled model.     

A theory of finite tensile deformation of double-network hydrogels

A continuum damage model was developed to describe the finite tensile deformation of tough double-network (DN) hydrogels synthesized by polymerization of a water-soluble monomer inside a highly crosslinked rigid polyelectrolyte network. Damage evolution in DN hydrogels was characterized by performing loading-unloading tensile tests and oscillatory shear rheometry on DN hydrogels synthesized from 3-sulfopropyl acrylate potassium salt (SAPS) and acrylamide (AAm). The model can explain all the mechanical features of finite tensile deformation of DN hydrogels, including idealized Mullins effect and permanent set observed after unloading, qualitatively and quantitatively. The constitutive equation can describe the finite elasto-plastic tensile behavior of DN hydrogels without resorting to a yield function. It was showed that tensile mechanics of DN hydrogels in the model is controlled by two material parameters which are related to the elastic moduli of first and second networks. In effect, the ratio of these two parameters is a dimensionless number that controls the behavior of material. The model can capture the stable branch of material response during neck propagation where engineering stress becomes constant. Consistent with experimental data, by increasing the elastic modulus of the second network the finite tensile behavior of the DN hydrogel changes from necking to strain hardening.     

移动加热器法生长碲锌镉晶体的组分输运与界面形貌研究

移动加热器法(THM)生长碲锌镉晶体时,界面稳定性对晶体生长的质量有很大影响。本文基于多物理场有限元仿真软件Comsol建立了THM生长碲锌镉晶体的数值模拟模型,讨论了Te边界层与组分过冷区之间的关系,对不同生长阶段的物理场、Te边界层与组分过冷区进行仿真研究,最后讨论了微重力对物理场分布的影响,并对比了微重力与正常重力下的生长界面形貌。模拟结果表明,Te边界层与组分过冷区的分布趋势是一致的,在不同生长阶段,流场中次生涡旋的位置会发生移动,从而导致生长界面的形貌随着生长的进行发生变化,同时微重力条件下形成的生长界面形貌最有利于单晶生长。因此,在晶体生长的中前期,对次生涡旋位置的控制和对组分过冷的削弱,是THM生长高质量晶体的有效方案。     

Nonlinear elastic constitutive relations of residually stressed composites with stiff curved fibres

Applied Mathematics and Mechanics - Mechanical models of residually stressed fibre-reinforced solids, which do not resist bending, have been developed in the literature. However, in some residually...     

A posteriori error analysis of a quadratic finite volume method for nonlinear elliptic problems

In this article, we construct and analyze a residual-based a posteriori error estimator for a quadratic finite volume method (FVM) for solving nonlinear elliptic partial differential equations with homogeneous Dirichlet boundary conditions. We shall prove that the a posteriori error estimator yields the global upper and local lower bounds for the norm error of the FVM. So that the a posteriori error estimator is equivalent to the true error in a certain sense. Numerical experiments are performed to illustrate the theoretical results.     

基于2,5-二溴对苯二甲酸配体的锌/钴配合物合成、晶体结构及性质

将有机物2,5-二溴对苯二甲酸(H₂L₁)和2,2′-联吡啶(L₂)作为双配体,使用溶剂热法和七水合硫酸锌(ZnSO₄·7H₂O)、六水合硝酸钴(Co(NO₃)₂·6H₂O)分别反应,得到配合物[Zn(L₁)(L₂)(H₂O)]_n(1)和配合物[Co(L₁)(L₂)(H₂O)]_n(2)。采用单晶X射线衍射、元素分析、红外光谱、紫外光谱、荧光光谱、热重分析等测试方法对这两种物质进行分析研究。单晶测试结果表明配合物1是单斜晶系,以Zn²⁺配位连接L^2-₁与L₂形成一维链状结构,各条链在分子间氢键和π…π共轭作用下有规律地堆叠形成三维网络结构。配合物2是三斜晶系,Co1离子和Co1ⁱ离子由H₂L₁上的羧酸氧原子O4和O4ⁱ连接,形成双齿螯合的配位结构单元,以Co²⁺配位连接 L^2-₁和L₂形成二维网格结构,各层在O—H…O分子间氢键和范德瓦耳斯力作用下有规律的堆叠形成三维网络结构。配合物1和2均含有芳香杂环、羧基杂环和氮杂环,具有良好的荧光性质和热稳定性,最大发射波长分别为345 nm和333 nm。     

An efficient edge based data structure for the compressible Reynolds-averaged Navier–Stokes equations on hybrid unstructured meshes

An efficient edge based data structure has been developed in order to implement an unstructured vertex based finite volume algorithm for the Reynolds-averaged Navier–Stokes equations on hybrid meshes. In the present approach, the data structure is tailored to meet the requirements of the vertex based algorithm by considering data access patterns and cache efficiency. The required data are packed and allocated in a way that they are close to each other in the physical memory. Therefore, the proposed data structure increases cache performance and improves computation time. As a result, the explicit flow solver indicates a significant speed up compared to other open-source solvers in terms of CPU time. A fully implicit version has also been implemented based on the PETSc library in order to improve the robustness of the algorithm. The resulting algebraic equations due to the compressible Navier–Stokes and the one equation Spalart–Allmaras turbulence equations are solved in a monolithic manner using the restricted additive Schwarz preconditioner combined with the FGMRES Krylov subspace algorithm. In order to further improve the computational accuracy, the multiscale metric based anisotropic mesh refinement library PyAMG is used for mesh adaptation. The numerical algorithm is validated for the classical benchmark problems such as the transonic turbulent flow around a supercritical RAE2822 airfoil and DLR-F6 wing-body-nacelle-pylon configuration. The efficiency of the data structure is demonstrated by achieving up to an order of magnitude speed up in CPU times.